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and Z bosons will be produced copiously at the LHC proton-proton collider. We 
study the parton distribution dependence of the total production cross sections and rapid- 
ity distributions, paying particular attention to the uncertainties arising from uncertainties 
in the parton distributions themselves. Variations in the gluon, the strong coupling, the 
sea quarks and the overall normalisation are shown to lead to small but non-negligible 
variations in the cross section predictions. Ultimately, therefore, the measurement of these 
cross sections will provide a powerful cross check on our knowledge of parton distributions 
and their evolution. 



1. Introduction 

A precise knowledge of parton distribution functions (pdf 's) is absolutely vital for reliable 
predictions for signal and background cross sections at the LHC Uncertainties can arise both 
from the starting distributions, obtained from a global fit to deep inelastic scattering and 
other data, and from DGLAP evolution to the higher Q 2 scales typical of LHC hard scattering 
processes. 

There are several reasons why it is very difficult to derive 'one sigma' errors on parton 
distributions of the form ±5/j. In the global fit there are complicated correlations between a 
particular pdf at different x values, and between the different pdf flavours. To give an example: 
the charm distribution is correlated with the gluon distribution, the gluon distribution at low x 
is correlated with the gluon at high x via the momentum sum rule, etc. Secondly, many of the 
uncertainties in the input data or fitting procedure are not 'true' errors in the statistical sense. 
To give an example: the uncertainty in the high-x gluon in the MRST fits derives from a 
subjective assessment of the impact of 'intrinsic on the prompt photon cross section. 

However, faced with the difficulties in trying to formulate global pdf errors, what one can do 
is make a detailed assessment of the pdf uncertainty for a particular cross section of interest, by 
determining which partons contribute and at what x and Q 2 values, and then systematically 
tracing back to the data sets that constrained the distributions in the global fit. Individual pdf 
sets can then be constructed to reflect the uncertainty in the particular partons determined by a 
particular data set (see Section 2). The purpose of this paper is to show how this can be done for 
the particular (test) case of W and Z production at the LHC. The advantage of choosing these 
cross sections is that the theoretical technology for calculating them is very robust. The total 
cross sections are known to NNLO in QCD perturbation theory and the input electroweak 
parameters (Mw,z, weak couplings, etc.) are known to high accuracy. The complete 0(a) 
electroweak radiative corrections are known, see for example Ref. [|J.Q Theoretically, therefore, 
the main uncertainty derives from the input pdfs and, to a lesser extent, as-f\ 

Knowledge of the uncertainty on the theoretical prediction allows one to assess the impact of 
an eventual experimental measurement. For example, with what precision would a measurement 
of, say, o~w have to be made to provide an additional constraint on the pdfs? We will not address 
any experimental issues here. Some studies of the likely experimental precision have already 
been performed, see for example Ref. with encouraging results. 

Before discussing the W and Z cross sections in detail, we begin with some general remarks 
about the range of x and Q 2 values that are probed by LHC cross sections. Figure |1] shows the 
values corresponding to the production of a heavy object (e.g. a W or Higgs boson, a tt pair, 
a multijet final state etc.) of mass M and rapidity y. We assume leading order kinematics, so 
that x = M exp(±y) / \fs and Q = M. As an example, a W boson (M = 80 GeV) produced 

1 Since the electroweak radiative corrections depend on the nature of, and experimental cuts on, the final 
state W and Z decay products H , we do not include them in our study. 
2 The two are of course correlated, see for example W. 
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at rapidity y = 3 corresponds to the annihilation of quarks with x = 0.00028 and 0.11, probed 
at Q 2 = 6400 GeV 2 . Notice that quarks with these x values are already 'measured' in deep 
inelastic scattering (at HERA and in fixed-target experiments respectively), but at much lower 
Q 2 , see Fig. § 

2. Parton sets used in this study 

In Ref. PJ we presented a set of pdfs, labelled MRST, that gave the best overall fit to the 
global set of data in that analysis. Associated with this set was a 'central' gluon distribution 
(at the starting scale of Qq — 1 GeV 2 ) and a 'central' value of as(M^) = 0.1175. This set 
resulted from a NLO analysis in the MS scheme and we can regard this particular set as a 
base from which we wish to explore the degree of latitude that relevant parameters can range 
over while still maintaining an acceptable description of all the data. Already two degrees of 
freedom were explored in Ref. [jl], (i) the variation of the average transverse momentum (kx) 
in prompt photon production, which resulted in a range of gluon distributions at large x, and 
(ii) the variation in the value of as(M§). These uncertainties were reflected in alternative sets 
of pdfs MRST(g |), MRST(g j), MRST(a 5 |T) an d MRST(a 5 [[) in which the extremes of an 
allowed range of values are taken and still yield acceptable fits. Again these sets were relevant 
to a NLO analysis in the MS scheme.0 

As we shall see below, the level of precision with which we should be able to make predictions 
for W and Z cross sections at the Tevatron and LHC is of the order of ±5%. It goes without 
saying that the numerical precision with which the pdfs are extracted and evolved should be 
much smaller than this. It is therefore worrying that a comparison of the aw values obtained 
using the default MRST partons and the recent (default) CTEQ5 partons |J reveals differences 
at the level of 5%, even though the data sets, theoretical assumptions and fitting procedures 
are largely the same in the two analyses. This difference has prompted us to re-examine the 
numerical accuracy of the MRST evolution code. As a result we discovered two small errors 
which affect the resulting pdfs at the ~ 1% level. Further details can be found in the Appendix. 
The basic five parton sets used in this analysis are not therefore the five sets presented in Ref. U, 
but rather their corrected versions, which we denote by MRST99, MRST99(g |), . . . etc.[] The 
W and Z Tevatron total cross section predictions using the MRST99 sets are some 2% larger 
than the corresponding MRST values presented in Ref. |]]|] The difference between the CTEQ5 
and MRST cross sections is only slightly reduced by using the corrected MRST99 partons. We 
will return to this below. 

Varying the gluon and as about their central values does not span the full range of uncer- 
tainties in the pdfs. Obviously, varying the normalisation of sets of input data to the analysis 



3 In ref. || we extended the analysis to provide equivalent sets of pdfs in the DIS scheme and also at LO. 
4 The FORTRAN code for all the MRST99 parton sets described in this paper can be obtained from 



frrttp : //durpdg . dur . ac . uk/HEPDATA/PDF , or by contacting W.J.Stirling@durham.ac.uk. 



5 In contrast, the LHC cross sections are essentially unchanged, due to a slight x dependence of the correction 
factor. 
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will result in a corresponding variation of the output pdfs and so this is an additional degree of 
uncertainty which we now explore. Since in this work we are primarily interested in uncertain- 
ties which most affect W and Z production at collider energies, i.e. uncertainties at small x, we 
focus mainly on the normalisations of the HERA structure function data. The quoted normal- 
isation uncertainties of F 2 from HI and ZEUS vary with Q 2 , but a mean value of 2.5% 0, || is 
appropriate. We then proceed to make two fits to all the data where we renormalise the HERA 
data first upwards by 1.025, and then downwards by 0.975, and allow the normalisations of the 
remaining DIS datasets to vary within their quoted uncertainties in order to get the best fit. 
We label these two pdf sets as MRST(g |) and MRST(g |). 

In Figs. and |3| we show comparisons of these two pdf sets with the default set at two 
values of Q 2 , a low Q 2 value typical of the Q 2 range of the fitted DIS data, and a high Q 2 
value typical of that probed in W, Z production at hadron colliders. One can see that when 
u and d increase to describe the larger F 2 , momentum conservation demands a decrease in the 
gluon distribution and vice versa. Because of the way charm and bottom pdfs are generated (by 
photon-gluon fusion), the heavy flavour pdfs simply 'mirror' the modified gluon. Increasing the 
normalisation also produces, at large x, a slightly increased d/u ratio (and vice versa) but this 
corresponds to a very small variation of F% '/F^. For example, at x — 0.675 and Q 2 = 35 GeV 2 , 
F^jFl = 0.447 ± 0.008. 

This variation of the normalisation effectively implies that the pdf sets MRST(g f ) and 
MRST(g I) give a measure of overall uncertainty on the u and d quarks. Next we try to get 
a corresponding measure of uncertainty on the strange quark distribution. The strange sea is 
assumed to be identical in shape to the non-strange sea at the starting scale Qq, but there 
is a strong flavour asymmetry due presumably to the larger mass of the strange quark. A 
suppression factor 2s/(u + d), at Qq — 1 GeV 2 of 0.5 is consistent with the di-muon production 
data of CCFR [[J. Fig. [| shows that a ±10% variation of this suppression is a realistic measure 
of the uncertainty in the CCFR estimate. Consequently we provide two more sets of pdfs 
corresponding to these upper and lower estimates for the strange sea suppression, labelled 
MRST(s T) and MRST(s |) 

The only parameter characterising the charm distribution in the MRST analysis is the charm 
mass, m c . In the updated default set m c = 1.43 GeV and we find a satisfactory description of 
F 2 C provided the variation of m c is within 0.15 GeV [ID|. So to complete our family of pdfs which 



explore individual degrees of uncertainty, we provide two sets MRST(c f) and MRST(c |) in 
which m c = 1.28 and 1.58 GeV respectively. A summary table of this 'MRST family' is given 
in Table 1. Note that an additional set MRST(d/tt f)) will be introduced and described in 
Section 4. 

3. W and Z total cross sections 

The total W and Z cross sections probe quark distributions in a fairly narrow x region 
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set 


a s (Ml) 

*J V Zj / 


comment 


MRST 


0.1175 


default set 


MRST(g t) 


0.1175 


larger gluon at large x 


MRST(# |) 


0.1175 


smaller gluon at large x 


MRST(a 5 tt) 


0.1225 


larger as 


MRST(a 5 ||) 


0.1125 


smaller as 


MRST(g t) 


0.1178 


larger (~ +2.5%) quarks at large x 


MRST(g |) 


0.1171 


smaller (~ —2.5%) quarks at large x 


MRST(s t) 


0.1175 


larger (~ +10%) strange quarks 


MRST(s |) 


0.1175 


smaller (~ —10%) strange quarks 


MRST(c t) 


0.1175 


larger charm quarks (m c = 1.28 GeV) 


MRST(c |) 


0.1175 


smaller charm quarks (m c = 1.58 GeV) 


MRST(d/u t) 


0.1175 


larger d/u ratio at large x 



Table 1: The various MRST99 pdfs used in the present study (see the Appendix). Note that 
the statements about the variations of the different quarks and gluons are relative to the default 
MRST99 set and apply at or just above the starting scale Ql = 1 GeV 2 . 

around x ~ Mv/Vs and at Q 2 ~ M£ ~ 10 4 GeV 2 .Q The inclusive rapidity distribut ions da/dyy 
select particular values of parton x (see Fig. |l|) but, at least in the case of W ± production, are 
presumably more difficult to measure. They will be discussed in the following section. 

We have already argued that the dominant theoretical uncertainty comes from parton dis- 
tributions, and so we fix all electroweak parameters at standard values [] and vary only the 
input pdfs (and as where appropriate). 

First, we consider the parton flavour decomposition of a w +, a w - and oz- Since the cross 
sections are dominated by the leading-order qq — > V processes, it is sufficient to consider the 
relative contributions of these to the total. Figures || and |6] show how the relative contributions 
of the various qq processes change with collider energy. We split the collider energy range at 
yfs = 4 TeV, and assume proton-antiproton and proton-proton collisions below and above this 
value. For production (Fig. H) we see that the ud W + and du — > W" contributions 
dominate at all collider energies. These are dominated by valence-valence and valence-sea scat- 
tering at pp and pp colliders respectively. The next largest contributions come from cs —>■ W + 
and sc — > W~ . Although these are sea-sea processes, they dominate the Cabibbo suppressed 
us — ► W + etc. valence-sea contributions. Perhaps the most significant result from Fig. |5] is 

6 For the subset of Ws and Zs produced with large transverse momentum the next-to-leading order qg 
scattering processes are also important, and as a result there is additional sensitivity to the gluon distribution 
in roughly the same x region - we will not discuss this further here. 

7 The Standard Model values B(W -> Iv) = 0.1084 and B(Z -> l + l~) = 0.03364 are used. The electroweak 
boson masses are taken to be My/ = 80.43 GeV and Mz = 91.1887 GeV. 
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that strange-charm scattering is a relatively large component of the total cross sections at LHC 
- approximately 17% for W + and 23% for W~ . This contrasts with their 5% contribution at 
the Tevatron. The remaining scattering processes contribute between 1% and 3% at LHC. 

The corresponding situation for Z° production is somewhat simpler, see Fig. |[ At LHC 
energies we have uu ~ dd ^> ss ^> cc, in line with the ordering of the parton distributions 
at the relevant x and Q 2 values.0 At the Tevatron, the ss and cc contributions are relatively 
smaller, and the uu contribution is correspondingly larger. 

We have already shown in Ref. |TJ that the Standard Model predictions for a(W) and 
cr(Z) at the Tevatron (y/s = 1.8 TeV) using the default MRST98 pdf set are in excellent 
agreement with the CDF and DO measurements. Here we repeat the comparison now including 
the extended set of MRST99 pdfs described in the previous section (and the Appendix). From 
Fig. |^, we would not expect any significant effect from varying the charm and strange partons, 
and indeed the W and Z cross sections calculated with the st,s|,c|,c| sets agree with the 
default MRST99 cross section to within less than one percent. On the other hand, the effect 
of changing the overall quark normalisation should be directly reflected in the cross sections. 
Figs. | shows the predictions of the first seven MRST99 sets (the default MRST99 set, together 
with six other X J sets (X = g,a s ,q) listed in Table 1) for a w and o z , together with the 
experimental data, all normalised to <r(MRST)|] There is an overall spread of approximately 
±3% about the default prediction,^ significantly smaller than the current experimental errors. 
The variations in the predictions are easily understood. At these (small) x values, the quark 
distributions increase with increasing /i 2 . The larger the as the faster the increase, and so 
°v( a s II) < crv(MRST) < cry(«5 tt)- The different gluon distributions also give rise to 
differences in the ay predictions. In this x ~ 0.05 region, the ordering of the gluon distributions 
is g | < g(MRST) < g I, see Fig. 2 of Ref. |l[ The larger the gluon the more rapid the DGLAP 
evolution, and so oy{g f) < ay(MRST) < oy{g |). The gluon variation is slightly smaller 
than the as variation. The largest variation (approximately ±3%) comes from the overall 
quark normalisation. Naively, one might have expected a variation twice as large, reflecting 
the variation in the q f j u and d partons near the starting scale, but the effect of DGLAP 
evolution to Q 2 ~ 10 4 GeV 2 is to suppress the differences. 

We may conclude from Fig. || that the net uncertainty in the &w,z predictions at the Tevatron 
is no more than about ±4% and is largely due to the normalisation uncertainty in the input u 
and d distributions, although DGLAP evolution dilutes this effect by about a factor of 2. 

8 Note that the weak neutral couplings are slightly different for u- and d-type quarks, and are taken into 
account in the curves shown in Fig. ^. 

9 Note that an important component of the experimental error is due to the luminosity measurement and 
uncertainty. The latter is quoted as ±3.6% for CDF and ±5.4% for DO 0. In addition, the value assumed 
for the total pp cross section is slightly different for the two experiments, and this may account in part for the 
systematically smaller DO cross sections displayed in Fig. ||. 

10 A measure of the scale dependence of these predictions is obtained by using instead fi = My/2 and 2My. 
The effect is shown as error bars on the default prediction and is evidently very small. 
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Figure [| shows the corresponding total cross sections at the LHC pp collider. We have now 
also included the predictions from the s j j,c j J, sets, but again we see that the effect on 
the total cross section is minimal. Even though the strange-charm contribution is relatively 
more important, DGLAP evolution has reduced the variation in the partons themselves. Inter- 
estingly the variations in the cross sections due to the gluon and the quark normalisation are 
approximately the same as at the Tevatron, but now the largest variation comes from as J. the 
reason for this is that DGLAP evolution is more rapid at the smaller x values probed at the 
LHC: schematically we have dq/dlogQ 2 ~ otsC where the (positive at small x) coefficient C 
increases as x decreases. Therefore the effect on the evolved q(x, Q 2 ) of varying as is magnified 
as x decreases. Bearing in mind that our as variation is somewhat conservative (the 'world 
average' error is now regularly quoted as ±0.004 or even less), it would appear from Fig. ||that 
±5% is a reasonable upper limit on the pdf uncertainty on aw,z at LHC. 

In Fig. | the (unequal) W + and W~ cross sections have been added together. It is interesting 
to consider what can be learned from the difference between these, in particular from the ratio 
Rzp = <Jw- / <J w+- Ignoring the contributions from strange-charm scattering (which contribute 
equally to aw+ and aw-), we have 

_ du d u 

R T «-==---= . (1) 
ud u d 

Notice that in the region of x we are considering, both fractions on the right-hand side are < 1. 
This raises the question of whether we can obtain any information on the ratio u/d, assuming 



that d/u is already well pinned down. Figure [10] shows the (proton-proton) collider energy 



dependence of the cross section ratio calculated using MRST (solid curve). For comparison we 
also show (dashed curve) the predictions for a modified version of MRST in which u is set equal 
to c?.0 Unfortunately the difference is small, and decreases with increasing yfs. The reason for 
this is that in the MRST fits, it is assumed that the asymmetry (d — u)/(d+ u) goes rapidly to 



zero as x — > 0. This trend is certainly suggested by the E866 Drell-Yan asymmetry data fl5 
However these data extend down only to x ~ 0.05, well above the x region probed by aw at 
LHC. In principle one could investigate the effect of 'pathological' behaviour of d — u as x — > 
on R T . 

At y/s = 14 TeV, the predictions for R T using all twelve MRST sets listed in Table 1 lie in 
the range 

R T = 0.731 ± 0.005 (2) 

with the central value corresponding to the MRST prediction. Evolution, normalisation etc. 
effects all tend to cancel in the ratio, leaving this as a 'gold-plated' prediction of the Standard 
Model. 

Before considering the W and Z rapidity distributions at the LHC, we consider the total 
cross section predictions using the latest CTEQ5 parton distributions |§. Since these are ob- 
tained in a similar way to the MRST partons, we would expect the (default, at least) predictions 



ii 



Formally, we replace the u and d pdfs appearing in the cross section expressions by their average. 
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to lie well within the MRST99±5% band. Surprisingly, this is not the case. Figures 11 and 12 
show the CTEQ5M and CTEQ5HQ 12 predictions compared to those of the various MRST99 
sets. The CTEQ5M predictions are about 4% (7%) larger at the Tevatron (LHC) than those 
of MRST99! The CTEQ5HQ predictions are even larger, lying well outside the ±5% band. 

The differences between the MRST99 and CTEQ5 cross section predictions are due to 
corresponding differences in the underlying quark distributions. This is illustrated in Fig. O, 



which shows the ratio of CTEQ5HQ and MRST99 partons at two representative Q 2 values, 10 
and 10 4 GeV 2 . The most striking difference is in the gluon at large x. The MRST and CTEQ5 
analyses put more emphasis in the global fit on the high (fixed target) prompt photon and 
(Tevatron collider) jet data respectively, and this is directly reflected in the resulting gluon pdfs. 
Since the sea quark pdfs are driven by g — > qq in the DGLAP evolution, the same behaviour 
is seen in the ratio of the strange and charm distributions in Fig. Of more relevance in 

the present context are the differences in the u and d quarks for x < 0.1 at Q 2 = 10 4 GeV 2 . 
One might expect that the long evolution length would result in essentially identical quark 
distributions^ but the CTEQ5HQ u and d quarks are between 2% (at x ~ 0.1) and 5% (at 
x ~ 0.0001) larger than those of MRST99 at small x. A possible reason for this discrepancy 
can be found in the numerical results presented in the Appendix. 

4. W and Z rapidity distributions 

The W + and W~ rapidity distributions are symmetric about yw = and have a distinctively 
different shape at large yw, see Fig. The bump in the W + distribution around yw = 3 

is caused by a very small x d quark scattering off a large x valence u quark. At large yw, 
therefore, the rapidity distributions probe small x (evolved) 'HERA quarks' scattering off large 
x (evolved) 'fixed target' quarks. For W + and W~ production, the latter are constrained by, 
for example, F 2 V and F^ structure functions respectively. 

The distributions themselves are presumably difficult to measure accurately in practice — it 
is the lepton rapidity distribution that is measured instead. We will not discuss the experimental 
considerations here, see for example [Q for more discussion. However it is worth pointing out 
that the ratio of rapidity distributions is sensitive to the ratio of d and u quarks. Thus 

dcr/dyw(W~) d(xi)u(x 2 ) d(xi) ufa) 



w ) 



da / dy w (W+) u(xi)d(x 2 ) u(xi) d(x 2 ) 



12r 



2 The 5M and 5HQ sets differ in their treatment of the heavy quark distributions ||, with the 5HQ distribu- 
tions more comparable to the MRST treatment Q. In order to maintain the same charm contribution to the 
F2 structure function at small x, the 5M and 5HQ gluon distributions are slightly different in this region, and 
this leads to corresponding differences in the W, Z cross sections, cf. the MRST g J variation. 

13 The smaller choice of m c = 1.3 GeV in the CTEQ5 analysis, compared to m c = 1.43 GeV for MRST99, 
also enhances the difference in the charm distributions. 

14 Note that at the lower Q 2 value in Fig. [l^ the impact of the different MRST and CTEQ5 starting param- 
ctcrisations and Q 2 cuts on the fitted DIS data can still be seen. At the higher Q 2 value these differences have 
evolved away and the ratios are more uniform. 

15 The curves are calculated using NLO corrections only - the NNLO corrections are not known at present. 
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« ^4 > x i = ex P(^) = ( 5 - 74 x 10 " 3 ) ex P(^) > ( 3 ) 

u{Xi) V s 

where for simplicity we have ignored the contributions involving strange and heavier quarks. 
Assuming that d/u — > 1 as x — > 0, the \axge-yw rapidity distribution ratio is therefore a 
direct measure of the large-x d/u ratio of partons, at Q 2 = M w . Interestingly, there has been 
considerable recent interest in this ratio, since it is also measured by charged current scattering 
cross sections at HERA: 

d 2 a/dxdQ 2 (e + p — > vX) 2 d(x,Q 2 ) 

d 2 a/dxdQ 2 (e-p vX) ~ (1 ~ V > u(x,Q 2 ) ^ 



Fig. [15| shows the ratio R^(yw) as a function of yw- Note that only the predictions of 
the MRST set are shown - the other ten sets of Table 1 give almost identical predictions. 
This is because the J variations quantified by these other sets largely cancel in the ratio. Also 
shown (dotted line) is the d/u ratio evaluated at x — (5.74 x 10~ 3 ) exp(yw) and Q = Mw- 
The convergence of the curves at large yw is confirmation of the derivation in (|3]). At central 
rapidities, the parton ratio falls below the W cross section ratio. This is due to the non-(u, d) 
contributions and u ^ d at the higher x values probed. For example, Fig. |14] shows the rapidity 
dependence of the strange-charm scattering contribution, which contributes equally to W + and 
W~ production. The W bosons from this scattering process are confined to the central 

region. 

As mentioned above, current information on d/u at large x comes from fixed target F^ 1 j 'F% v 
measurements and the lepton asymmetry in pp — > W ± + X, see for example In the MRST 
fit, NMC n/p data are used to constrain the large-x <i-quark pdf in this way. When the 
corresponding predictions for &cc( e+ p) are compared with the ZEUS data |16| , there is some 
evidence of a slight excess of data over theory in the largest x (= 0.42) bin. Could this imply 
that the d/u ratio is being underestimated in the standard global fits? Any attempt to increase 
d/u at large x in the global fit leads to a direct conflict with the n/p data. However, Bodek and 
Yang |17| have recently stressed a point first made by Thomas and Melnitchouk (Ref. |18| and 
references therein) that the n/p data should be corrected for nuclear binding effects which, at 
large x, lead to a larger d/u ratio, in 'better' agreement with the ZEUS data. This is an issue 
that has yet to be resolved, and therefore independent information from the W ± rapidity ratio 
at LHC would be very valuable. 



To illustrate the likely size of the effect, we show in Fig. [15] the predictions of a modified 
version of MRST in which the d/u ratio is deliberately enhanced at large x such as to mimic 
the size of effect proposed in Ref. [T7|]. The difference is only significant for y w > 4, which 
suggests that discrimination will not be easy at the LHC. 

Since it is the charged lepton rapidity distribution that is measured in practice, we can see 
to what extent the differences between W + and W~ production survive the leptonic decay 
process. Fig. [16] shows the ratio RL of the charged lepton rapidity distributions.^ A transverse 

16 The lepton distributions arc computed at leading order with the full V — A decay structure taken into 
account. The ratio is very insensitive to higher-order effects. 



s 



momentum cut of p l T > 25 GeV/c is imposed. The solid line is again the MRST prediction, and 



the dashed line is the prediction of the enhanced d/u set, as in Fig. [15]. Notice the qualitatively 
different behaviour of the ratio at large rapidity. This is due to the V — A structure of the 
qq' — > tvi subprocess. Thus a W + produced with large positive yw in ud — > W + , decays 
preferentially to a e + with smaller rapidity ye, while a W~ produced with large positive yw in 
du — > W~, decays preferentially to a e~ with larger rapidity ye- The ratio i?l therefore increases 
as ye —>■ y max - This is in direct analogy to the lepton asymmetry in pp — > WX changing sign 
and becoming negative at large yg. The effect of enhancing the d/u ratio at large x is still just 
visible. 

In the same spirit we can consider the ratio of and Z° rapidity distributions. The 
analogue of (|3|) is now 

R / x da/dy(Z°) k u ujx'Jujx^) + K d d(x' l )d(x' 2 ) 

Z/W[V) ~ da/dy(W+) + da/dy(W-) ~ \ V J? {u( Xl )d(x 2 ) + d{x 1 )u(x 2 )} ' 1 ' 

Here the u and d contributions to the numerator are weighted by the appropriate neutral current 
coupling factors: 

K q = cos" 2 9 W (v 2 q + a 2 ,) w 0.37(0.48) for g = u(d) . (6) 

Because M\y ^ Mz, the pdfs in the numerator and the denominator are no longer evaluated 
at exactly the same values of x. This means that there is no simple analytic expression for the 
ratio in the y — > y max limit. The prediction for Rz/w{y) calculated using the MRST partons 
is shown in Fig. [I7[ The prediction using the MRST^/ii-enhanced) set is essentially the same 
— there is little sensitivity to the d/u ratio. 



5. Conclusions 

We have made a detailed study of the W and Z boson total production cross sections and 
rapidity distributions at the LHC proton-proton collider. We have focused on the dependence 
of the cross sections on the parton distributions — the largest source of theoretical uncertainty. 
We have shown that the two important effects are the DGLAP evolution uncertainty, due to the 
uncertainty in as and to a lesser extent in the gluon, and the normalisation uncertainty on the 
DIS structure function data fitted in this analysis, which is directly reflected in a corresponding 
uncertainty in the parton distributions. Taking everything together, we conservatively estimate 
an overall pdf uncertainty of approximately ±5% on the total cross sections at the LHC, only 
slightly larger than the uncertainty on the corresponding Tevatron cross sections. 

We have also studied the rapidity distributions. If these can be accurately measured, they 
allow the quark distributions to be probed over a wide range of x. We showed, in particular, that 
the difference in the W~ and W + distributions at large rapidity directly reflects the difference 
in the d and u distributions as x — *■ 1. A measurement of the d/u ratio in this way would be 
complementary to the information obtained from charged current e ± p cross sections at HERA. 
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Finally, while comparing our total cross section predictions with those based on the latest 
CTEQ5 distributions we noticed differences that we believe reflect shortcomings in the accuracy 
of the evolution codes. We have corrected an error (at NLO) in the standard MRS evolution 
code that brings it into exact agreement with the "standard evolution" benchmark program. 
The CTEQ5 W and Z cross sections are still systematically larger than those using the corrected 
MRST partons, an effect we believe is due to a too rapid CTEQ sea quark evolution. 
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Appendix 



The parton sets used in this paper, MRST99, differ slightly from the original MRST [Ij] 
distributions. The main difference is due to the correction in the computer code of a factor 
Cp —>■ Ca in one of the terms appearing in the NLO P gg splitting function. The effect is to make 
the gluon, and hence the sea, evolve a little more rapidly at small x. We have also included a 
missing NLO P cq contribution which affects the charm evolution only. 

It is informative to explain the present situation concerning evolution codes. The long 
evolution length to Tevatron and LHC energies, and the expected precision of the predictions 
of the weak boson production cross sections, can expose small defects in evolution codes. The 
unexpectedly large difference that we found between the cross sections predicted using MRST98 
and CTEQ5 partons means that all codes should be carefully checked. The difference appears 
to be due to the more rapid evolution of the sea quarks of CTEQ5. A similar effect was found 
in a comparison of various NLO evolution codes made during the 1996 HERA workshop fl9 



Neither CTEQ, nor MRS, agreed completely with the "standard code" for a test evolution 
from Q 2 = 4 to 100 GeV 2 . In the interval 10" 4 < x < 10 -1 , the CTEQ sea was some 4% high, 
and the MRS sea was 1% low and the gluon about 1.5% low. After including the corrections 
noted above, the MRST99 code is found to agree precisely with the "standard" evolution, 
and consistent with exact momentum conservation as a function of Q 2 , as it should be. The 
momentum fractions carried by the various partons as Q 2 increases are shown in Table 1. Small 
violations of overall momentum conservation can be seen for MRST98 and the CTEQ parton 
sets (CTEQ4M and CTEQ4HQ are qualitatively similar). 

Since the MRST99 partons come from a new global fit to the data, a comparison with 
the MRST98 partons Q does not simply reflect the small correction to the evolution code. 
In particular, ZEUS 1995 SVX data |20| is included in the MRST99 analysis. Moreover we 
have increased the mass of the charm quark m c from 1.35 to 1.43 GeV so as to maintain the 
description of the F 2 (charm) data. 
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CTEQ5M 



u r 



2u x 



2d, 



2s 



2c 



2b 



9 



total 



Q 2 (GeV' 



= 5 
10 
50 
100 
500 
1000 
5000 
10000 



0.2843 
0.2667 
0.2374 
0.2278 
0.2096 
0.2031 
0.1902 
0.1855 



0.1143 
0.1072 
0.0955 
0.0916 
0.0843 
0.0817 
0.0765 
0.0746 



0.0559 
0.0586 
0.0631 
0.0646 
0.0674 
0.0684 
0.0703 
0.0710 



0.0730 
0.0746 
0.0774 
0.0783 
0.0800 
0.0806 
0.0818 
0.0822 



0.0374 
0.0412 
0.0477 
0.0498 
0.0538 
0.0552 
0.0579 
0.0589 



0.0120 
0.0174 
0.0265 
0.0294 
0.0350 
0.0370 
0.0409 
0.0424 



0.0000 
0.0000 
0.0059 
0.0097 
0.0169 
0.0194 
0.0245 
0.0263 



0.4293 
0.4424 
0.4576 
0.4607 
0.4667 
0.4688 
0.4727 
0.4742 



1.0063 
1.0081 
1.0110 
1.0120 
1.0137 
1.0142 
1.0149 
1.0150 



CTEQ5HQ 



d„ 



2u, 



24 



2s 



2c 



2b 



a 



total 



Q 2 (GeV- 



! ) = 5 
10 
50 
100 
500 
1000 
5000 
10000 



0.2858 
0.2681 
0.2386 
0.2290 
0.2107 
0.2042 
0.1912 
0.1864 



0.1160 
0.1088 
0.0968 
0.0929 
0.0855 
0.0828 
0.0776 
0.0756 



0.0579 
0.0604 
0.0646 
0.0660 
0.0686 
0.0695 
0.0713 
0.0719 



0.0749 
0.0763 
0.0787 
0.0796 
0.0811 
0.0816 
0.0826 
0.0830 



0.0379 
0.0416 
0.0479 
0.0499 
0.0538 
0.0552 
0.0579 
0.0589 



0.0118 
0.0171 
0.0260 
0.0290 
0.0345 
0.0365 
0.0404 
0.0418 



0.0000 
0.0000 
0.0058 
0.0096 
0.0167 
0.0192 
0.0242 
0.0260 



0.4215 
0.4354 
0.4519 
0.4554 
0.4621 
0.4645 
0.4689 
0.4705 



1.0058 
1.0075 
1.0104 
1.0114 
1.0130 
1.0135 
1.0142 
1.0143 



MRST98 



2u, 



2d, 



2s 



2c 



2b 



9 



total 



Q 2 (GeV* 



: ) = 5 
10 
50 
100 
500 
1000 
5000 
10000 



0.2799 
0.2627 
0.2341 
0.2240 
0.2061 
0.1994 
0.1868 
0.1818 



0.1160 
0.1089 
0.0970 
0.0928 
0.0854 
0.0826 
0.0774 
0.0753 



0.0599 
0.0614 
0.0641 
0.0652 
0.0670 
0.0677 
0.0690 
0.0695 



0.0759 
0.0764 
0.0775 
0.0779 
0.0787 
0.0790 
0.0796 
0.0799 



0.0414 
0.0441 
0.0487 
0.0504 
0.0533 
0.0545 
0.0567 
0.0575 



0.0087 
0.0130 
0.0205 
0.0232 
0.0281 
0.0300 
0.0335 
0.0349 



0.0000 
0.0000 
0.0062 
0.0096 
0.0156 
0.0179 
0.0222 
0.0239 



0.4108 
0.4245 
0.4406 
0.4452 
0.4524 
0.4551 
0.4594 
0.4609 



0.9926 
0.9910 
0.9887 
0.9882 
0.9866 
0.9862 
0.9846 
0.9838 



MRST99 



d v 



2u, 



2d, 



2s 



2c 



2b 



9 



total 



Q 2 (GeV* 



: ) = 5 
10 
50 
100 
500 
1000 
5000 
10000 



0.2810 
0.2640 
0.2359 
0.2259 
0.2083 
0.2016 
0.1891 
0.1841 



0.1136 
0.1067 
0.0953 
0.0913 
0.0842 
0.0815 
0.0764 
0.0744 



0.0600 
0.0616 
0.0646 
0.0657 
0.0676 
0.0684 
0.0697 
0.0703 



0.0760 
0.0767 
0.0780 
0.0786 
0.0795 
0.0799 
0.0805 
0.0808 



0.0416 
0.0444 
0.0492 
0.0509 
0.0540 
0.0552 
0.0574 
0.0583 



0.0073 
0.0122 
0.0205 
0.0235 
0.0287 
0.0308 
0.0345 
0.0360 



0.0000 
0.0000 
0.0066 
0.0102 
0.0166 
0.0190 
0.0235 
0.0253 



0.4197 
0.4335 
0.4482 
0.4521 
0.4584 
0.4608 
0.4648 
0.4662 



0.9992 
0.9991 
0.9983 
0.9982 
0.9973 
0.9971 
0.9960 
0.9953 



Table 2: Momentum fractions carried by quarks and gluon for the CTEQ5M ||, CTEQ5HQ ||, 
MRST98 U and the new MRST99 parton sets, for Q 2 between 5 and 10 4 GeV 2 . The fractions 
are obtained by numerically computing the integrals J^ o xfi(x,Q 2 )dx where xq = 1CT 5 . The 
missing interval, < x < 1CT 5 , accounts for the small departure of the MRST99 total from 1. 
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Figure 1: Values of x and Q 2 probed in the production of an object of mass M and rapidity y 
at the LHC, = 14 TeV. 
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Figure 2: Comparison of the MRST(g |) and default MRST pdfs at Q 2 = 10 and 10 4 GeV 2 
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Figure 3: Comparison of the MRST(g |) and default MRST pdfs at Q 2 = 10 and 10 4 GeV 2 
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Figure 4: Comparison of the strange sea from the NLO analysis of CCFR |§ dimuon produc- 
tion with the default MRST (solid line) and a variation of ±10% in the proportion (strange 
sea)/(total sea) at Qq = 1 GeV 2 representing the input to the pdf sets MRST(s f) and 
MRST(s I). 
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flavour decomposition of W cross sections 
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Figure 5: Parton decomposition of the W + (solid line) and W" (dashed line) total cross sections 
in pp and pp collisions. Individual contributions are shown as a percentage of the total cross 
section in each case. In pp collisions the decomposition is the same for W + and W~ . 
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flavour decomposition of Z cross sections 
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Figure 6: Parton decomposition of the Z° total cross sections in pp and pp collisions. Individual 
contributions are shown as a percentage of the total cross section. 
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Figure 7: Predictions for the total W, Z production cross sections times leptonic branching 
ratios in pp and pp collisions, as a function of the collider energy ^/s. The default MRST99 



partons are used. Experimental measurements from UAl [IT], UA2 |12 |, CDF [|13j and DO [|14 
are also shown. 
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W and Z Cross Sections: Tevatron 
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Figure 8: Predictions for the W and Z production cross sections times leptonic branching ratios 
in pp collisions at 1.8 TeV using the various MRST parton sets discussed in the text. The error 
bars on the default MRST prediction correspond to a scale variation of \i = My/2 — > 2My, 
V = W, Z. Experimental measurements from CDF and DO Jl4| are shown. 
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W and Z Cross Sections: LHC 
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Figure 9: Predictions for the W and Z production cross sections times leptonic branching ratios 
in pp collisions at 14 TeV using the various MRST parton sets discussed in the text. 
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Figure 10: Prediction for the ratio R T of W~ and W + total cross sections in pp collisions, as 
a function of the collider energy y/s. The default MRST99 partons are used. For pp collisions 
the ratio is I. Also shown (dashed line) is the prediction obtained by setting u = d in the quark 
sea. 
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W and Z Cross Sections: Tevatron 
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Figure 11: Predictions for the W and Z production cross sections times leptonic branching 
ratios in pp collisions at 1.8 TeV using the MRST99, CTEQ5M and CTEQ5HQ |J parton sets. 
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W and Z Cross Sections: LHC 
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Figure 12: Predictions for the W and Z production cross sections times leptonic branching 
ratios in pp collisions at 14 TeV using the MRST99, CTEQ5M and CTEQ5HQ @ parton sets. 
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Figure 13: Ratio of the partons of the CTEQ5HQ || set to those of the default MRST99 set 
at Q 2 = 10 and 10 4 GeV 2 . 
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Figure 14: Rapidity distributions of the W and W + bosons at the LHC. Also shown (dashed 
line) is the (common) charm-strange scattering contribution. 
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Figure 15: Ratio of the rapidity distributions of the W~ and W + bosons at the LHC. The 
solid line corresponds to the default MRST partons and the dashed line to a modified set (see 
Table 1) with an enhanced d/u ratio at large x. Also shown (dotted line) is the ratio of d to u 
pdfs evaluated at x = M w /yfs exp(y w ) and Q 2 



M 2 W . 
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Figure 16: Ratio of the charged lepton rapidity distributions from W~ and W + production and 
decay at the LHC. The solid line corresponds to the default MRST partons and the dashed line 
to a modified set (see Table 1) with an enhanced d/u ratio at large x. 
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Figure 17: Ratio of the rapidity distributions of the Z° and W ± bosons at the LHC, calculated 
using the default MRST partons. 
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